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The  Location,  by  Means  of  a  Single  Camera,  of 
A  Projectile  with  a  Known  Line  of  Fire 


Agree  that:  the  line  of  fire  is  tne  intersection  of 
the  vertical  trajectory  plane  with  the  horizontal  plane 
through  the  camera  (that  is  the  center  of  the  camera  lens) : 
the  line  of  eight  is  the  intersection  of  the  latter  plane 
with  the  vertical  plane  through  the  optical  axis;  s  is  the 
ingle  of  inclination  of  the  optical  axis;  cp  so  measures  the 
angle  formed  by  the  line  of  sight  and  the  perpendicular 
from  the  camera  to  the  line  of  fire  that  (p  Is  positive  in 
Fig*  £  and  negative  in  Fig.  II. 
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Thus  9  might  be  defined  as  the  amount  of  counterclockwise 
(picture  the  clock  face  up  and  directly  beneath  the  camera) 
turning  necessary  to  bring  the  line  of  sight  perpendicular 
to  the  line  of  fire,  it  being  understood  that  clockwise 
turning  is  accorded  a  negative  value.  Also  agree  that:  D 
is  the  distance  from  the  camera  to  tne  point  at  which  the 
line  of  sight  intersects  the  line  of  fire,  p  is  the  distance 
from  the  camera  to  the  point  0,  at  which  the  optical  axis 
pierces  the  trajectory  plane. 

Imagine  that  the  origin  of  3  space  is  so  placed  at  0 
that*  the  positive  first  axis  (say  k^)  is  horizontal  and 

passes  directly  above  the  camera;  the  second  (say  k2)  axis 
is  horizontal,  perpendicular  to  the  axis,  and  points  clock¬ 
wise  about  the  camera  (that  is  to  the  camera  man’s  right);  the 
third  (say  kj)  axis  points  vertically  upward.  Let  C  be  the  p 

point  at  which  we  find  the  center  of  trie  camera  lens  and  adopt 
the  general  convention  wnereby  c  3  (C^,  C2,  C-).  i*et  the  ^ 

axis  agrge  with  the  k~  axis  and  let  the  hi  axis  so  pass 
through  0  in  an  upward iy  direction  tnat  ^  is  perpendicular  to 
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both  the  optical  and  the  ^  axis*  Let  the  x*  axis  pass 

horizontally  through  the  origin  in  the  trajeotory  plane  and 
point  clockwise  about  the  camera,  and  let  the  y*  axis  agree 
with  the  k^  axis. 

It  should'  be  noted  that  a  ray  through  C  with  standard 
coordinates  (  5,  r) )  pierces  the  plane  determined  by  and 

^  in  a  point  whose  (^i7^)  coordinates  are  (p£,  prj). 

Let  u^1),  u^,  v^,  v<2>  be  the  uni*  points  on  the 
respective  positive  x*,  y*  axes,  ana  let 

T  *  v(l)  ^v(2). 


A  check  reveals: 

U(1)  =  (0,  1,  0), 

«W  *  (sin  e,  0,  cos  «), 

t 

v<l>  =  (-sin  <p,  cos  <p,  0),  i 

v(2)  =  (0,  0,  1), 

f  *  (cos  qp,  sin  <p,  0), 

C  =  (pcos  a,0,  -  psln  a). 

_j3ince  the  trajectory  plane  passes  thresh  v^, 

and  0,  it  i6  eviden  that  the  line  through  0  and  t  is  perpen¬ 
dicular  to  the  trajectory  plane.  Accordingly  the  trajectory 
plane  consists  of  those  phAnt*  K  for  which 
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Letting  P  be  any  point  and  Q  be  such  a  point  in  the 
trajectory  plane  that  Q  ia  g^Uinei&r  with  C  and  P,  choose  X 

so  that  . 


Clearly 


(1  -  X)C  V-XP  «  Q.  ‘"-V 

♦  r  ,  ’  *  '  •'  1  ■  %  '  l  ;»  •' 

(1  *  X  £ ♦  Xt  •  P  ■  -T  * 


K  .  I-S  ■ 

t  •  C  -  T  •  P 


T  -  C  -T •  p 


Q 


V 
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Suppose  nowt  that  P  and  Q  are  related  as  above j  that  I. 
the  photographic  image  of  Q,  has  plate  measurements  (p^,  p,j  j 

and  that  standard  coordinates  ( rj)  of  the  ray  CQ  have  been 
introduced  in  such  a  way  that 

i  P  *  *  prju^ 

I 

& 

f  .  5  *  A10  +  A11  P1  *  A12P2 

^  *  a20  *  a21  p2  +  a22  p2. 

fote  that  P  lies  in  the  plane  and  that  (p£,  p rj)  are 

the  coordinates  of  P  with  respect  to  the  r^)  axes. 

Let  (x*,  y*)  be  the  coordinates  of  0  with  respect  to  the 
horizontal  trajectory  x*  axis  and  the  vertical  y*  axis  and 
verify  thati 

X*  *  *  Q  y*  a  -  Qj 

X,  .  (T-  CHP-  T(lh  -  fl-PHC- 

y.  ,  ft  •  c) (p »  t <2b  -  ft.  phc-  v(2h 

*  T  »  C  T  •  P  * 

P  •  V^  a  pfc(u^  -  V^)  ♦  P  jjU^.  V<x> 

«  p£  cos  f  -  p  ijsin  cp  sin  t, 

p  .  yV2)  b  *  V(2>)  ♦  PT](u^2^'  v(2)) 


f 


t  •  p 


*  P  7}  COS  t, 

a  p  5  (U^*  T  )  4  p  rj(uv2)  •  t  ) 


*  p£  sin  (p  ♦  p  »]  cos  9  sin  c. 


T  *  C 


p  cos  9  cos  t, 


k 


* 


C  '  vtt)  =-p  sin  9  cos  *, 
C  .  v( 2 )  s  -  P  sin  c; 


(t  ♦  C  -  t  •  p)x*  *  (p  cos  9  cos  s)  |j»  5  cos  9  -  p  rj  sin  <p  sin  «J 
+  (p  sin  9  cos  a)  jj>  £,  sin  9  ♦  pn  cos  9  sin  s] 


=  &>  cos  a, 


(t  •  c  -  t  *  p)y*  3  (p  cos  9  cos  c) p  n  cos  c 

♦  (p  sin  .){p*  «ln  9  +  p  t\ cos  9  sin  ej 


2  2 
*  sin  9  sin  a  ♦  qp  c°*  9$ 


(T.  C  -  t*  P)  *  p  cos  9  cos  a  -  5P  sin  9  -  7$>  cos  9  sin  i) 


x*-  * 


U  co?„i. 


cos  9  cos  a  -  £  sin  9  -  J)  cos  9  sin  « 


y*  *  „S  A,slJa. <  *  nfl.Sfll.il _ 

cos  9  cos  •  -  5  sin  9  -  7)  cos  9  sin  e 


Now  let  F  be  the  intersection  of  the  horizontal  line  of 
tight  and  the  line  of  fire.  Let  the  x  axis  lie  along  the  line 
of  fire  and  point  clockwise  about  the  camera*  X**t'  the  y  axis 
pass  through  F  and  point  vertically  upward,  let  (x,  y)  be  the 
coordinates  of  Q  with  respect  to  these  a kes,  and  let  D  be  the 
dU  *nce  from  C  to  F.  It  follows  that 


D  *  p  cos  a,  y  3  y*  ♦  p  sin  a>,  x  *  x*. 
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x  s  cos  (p  cos  t  -  t,  sin  <p  -  rj  cos  cp  sin  c  * 

y  *  P-OO*  SL  %  +  nD  cos  9  cog  < 

*  cos  9  cos  e  -  5  sin  p  -  rj  cos  9  sin  &  * 

x  « _ LC_§S£_JL _ 

cos  s  -  £  tan  9  -  Tjsin  e  * 

v  « _ -jJ>.  * _  , 

*  cos  e  -  Ttan  <p  -  rjsin  e  ' 

*  “  A10  *  A11  PX  +  a12  p2' 

n  *  A20  *  A21  P1  +  A22  p2‘ 

■Hi®  last  two  equations,  which  are  fitted  as  a  rule  by 
least  squares,  formulate  the  passage  from  the  plate  measure¬ 
ments  (p^,  p2)  to  the  standard  coordinates  (5  ,rt  ). 

Some  computational  suggestions  follow.  1 


Agree: 

that  Matrix 

I  is 

cos  a 

-  tan 

<P 

-  sin  c 

0 

sec 

<P 

0 

sin  c 

0 

cos  e 

that  Matrix 

II  is 

1 

0 

0 

A10 

A11 

A12 

A2Q 

A21 

A22 
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that  Matrix  III  is  the  classical  row  into  column  result  of 
multiplying  Matrix  I  into  Matrix  II. 

Let 


be  Matrix  III, 


Note  that 


Roo 

R01 

R02 

R10 

R11 

R12 

R20 

and  let 

R21 

R22 

X  »  X 

A  P  » 

y  =  X 

x  *  DX, 

y  »  DY, 

and  that  the  computation  of  (X,  Y)  from  (p,,  p?)  may  be 
effected  in  the  following  way.  1 


Col¬ 

umn 

Entry 

Heading 

Foot¬ 

ing 

Checks 

1 

Frame  Numbers 

Frame  number 

2 

First  Plate 
Coordinates 

Plate  Measure¬ 
ment,  p2 

23 

3 

Second  Plate 
Coordinates 

Plate-  Measure¬ 
ment,  p2 

23 

4 

R00*R01P1+R02P2 

A 

24 

^snR00+R01?  2<fR02^ 

5 

(R10+R11p1*R12p2)/4 

X 

24-5 

24  •  4*nB10+R1:L  22^Rl2  23 

a 

(R20*R21 P1+R22P2^  ^ 

I 

24-3 

24  •  3=nR2Q+P.,-.  ^2^^^22  ^ 

N.B.  n  =  the  number  of  entries  in  coluan  1. 

V.B.  *12”  means  "sum  column  2";  "24. 5"  means  "sum  the  products 
pf  corresponding  entries  of  column  4  and  coluan  5." 


A  good  chock  on  the  matrix  multiplication  of  Matrix  I  1 
into  Matrix  II  may  be  based  on  the  inner  product  -  sum  formula 


2 

£  Ci 

J*0  * 


(2,2) 

1  *  0 

(M)*(o,o) 


where  r^  is  the  i^  row  of  Matrix  I,  is  the  J ^column  of 

Matrix  II,  and  R..  is,  as  agreed,  the  (±, J ) ^  element  of 
Matrix  III. 

If  the  above  reduction. is  used  on  the  present  (Spring  194 
data  (viz.  lake  473)  obtained  by  the  Photographic  Measurements 
8ection,  then  the  following  rules  should  be  followed: 

Use  5  decimal  trigometric  tables  and  make  all  entries  in 
Matrix  I  to  5  decimals,  as  in 


.99347 


.02892 

1.03020 


-.11411 


.11752 


.98366 


Make  entries  in  Matrix  II  to  5  decimals  in  column  1  and  to  7 
decimals  in  columns  2  and  3,. as  in 


.00251 


.0011949  .0000129 


.00065 


.0000023  .0011941 


Make  each  entry  in  Matrix  III  to  8  decimals,  as  in 
1,99345801  .00003429  -.00013589  || 


.00221493  .00123099  .00001329 


.11815938  .000o0226  .00117459 


frploy  the  Inner  product-sum  formula  to  obtain  a  7  decimal 
check.  Before  using  Matrix  III  round  each  of  its  entries  to  5 
decimals  in  column  1  and  to  7  decimals  in  columns  2  and  3, 

as  in 
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1  .99346 

.0000343 

-.0001359 

.00221 

.0012310 

.0000133 

.11816 

.0000023 

.0011746 

Matrix  III  may  now  be  used  in  obtaining  (X,  Y)  from  (p^,  p2) 

as  indicated  above.  In  this  connection  write  d,  X ,  and 
Y  to  5  decimals. 

Since  the  standard  error  to  be  feared  in  a  determination 
of  either  x  or  y  is  usually  about  .15  feet,  the  rounding  errors 
introduced  into  X  and  Y  should  be  negligible  as  long  as  D 
is  less  than  10,000  feet. 

Although  the  above  matrices  could  be  so  altered  that  their 
use  would  lead  directly  to  x  and  y  instead  of  X  and  Y  it  does 
not  appear  desirable  to  make  this  change.  Firstly,  frame  by 
frame  determination  of  (x,y)  is  seldom  needed.  In  fact  usually 
only  the  center  value  of  (x,y)  ia  of  Interest.  Secondly  it  is 
clear  that  any  linear  function  of  the  x's  and  y's  may  be  evalu¬ 
ated  by  applying  it  to  the  X‘s  and  Y's  and  multiplying  the 
result  by  the  scale  factor  D.  Thirdly  the  preheat  matrix 
arrangements  are  more  uniform  and  more  easily  remembered  than 
those  met  within  the  proposed  alternative. 

If,  as  is  frequently  the  case,  c  and  Matrix  II  remain 
constant  for  a  number  of  rounds  it  is  desirable  from  some 
points  of  view  to  arrange  the  computation  differently.  In 
Matrix  I  replace  "tan  f"  and  "sec  f"  by  "0"  to  obtain  Matrix  I_. 
Multiply  Matrix  I  into  Matrix  II  to  obtain  Matrix  III  and 


let  Bq,  B^,  B2  be  the 
■  (0,  0,  0)  recall 
and  observe 

three  rows  of  Matrix  III,..  Check  that 
the  distrlbutivety  of  matrix  multiplication 

(*W  R0l*  *  Bo 

♦  (-tan  f ) (A10*  ^12^ 

(R10»  Rll*  R12^  "  B1 

♦  (lee  f) (A^q,  l^l*  *12^ 

• 

(R20»  R21*  R22^  *  B2 

These  formulas  make  clear  the  simplicity  with  which  a 
calculating  machine  may  be  used  to  obtain  «ny  given  positional 
entry  in  Matrix  III  for  each  of  a  number  of  ronnds.  Simple 
sum  checks  suggest  the aselte*. 


* 


Xfcehould  be  noted,  however,  that  efficient  use  of  this 
method  depends  on  obtaining  Matrix  HI  for  each  of  a  number 
of  rounds  before  Matrix  lit  is  used  for  any  one  round. 


If  frame  by  frame  determination  of  x  and  y 
a  new  matrix  say  be  to  formed  that  its  rows  are 


are  desired, 
respectively 


Bq  ♦  (-  tan  $)(A^q»  ^il*  A^) 


(D  sec  tp)  A^i ,  A^) 


From  this  new  matrix  x  and  y  may  be  arrived  at  in  the  same 
manner  as  X  and  X  were  reached  via  Matrix  III. 


f\  .p.'WVr***. 
A'.  P. Morse 
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